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I.  SUMMARY  OF  RESEARCH 

Under  a  previous  contract  (N00019-77-C-0299)  an  aircraft  simulation 
model  was  developed  for  the  analysis  of  near  field  pattern  effects 
due  to  an  airborne  antenna.  The  solution  was  based  on  the  geometrical 
theory  of  diffraction  (GTO).  That  model  treated  the  fuselage  as  a 
composite  elliptic  cylinder  and  the  wings  and  stabilizers  as  finite 
flat  plates.  The  near  field  formulation  developed  has  been  used  to 
calculate  the  radiation  patterns  in  and  close  to  the  roll  and  elevation 
planes  of  the  aircraft.  The  numerical  results  obtained  were  in  excellent 
agreement  with  the  measured  patterns1. 

However  as  the  previous  solution  was  used  to  analyze  various 
private,  commercial  and  military  aircraft,  it  became  apparent  that 
our  simple  aircraft  model  was  unable  to  simulate  the  arbitrary  sub¬ 
structures  of  the  aircraft,  such  as  T-tails,  jet  engine-housing,  and 
numerous  structures  being  attached  to  the  basic  fuselage  and  wings. 

For  these  reasons,  the  previous  solution  was  extended  under  the  present 

p 

contract  (N00019-78-C-0524)  in  terms  of  a  multiple  plate  model  for 

the  aircraft  substructures.  This  multiple  plate  approach  allows  one 

to  model  the  various  substructures  by  boxing-out  the  structure  using 

a  set  of  finite  flat  plates.  The  finite  flat  plates  are  chosen  to 

be  the  basic  building  blocks  to  simulate  these  substructures  in  that 

they  are  easy  to  input  in  the  computer  code  and  are  efficient 

to  analyze.  The  basic  concept  of  the  multiple  plate  model  is  described 

here  and  one  specific  example  is  used  to  demonstrate  its  validity. 

More  numerical  results  with  accompanying  measurements  can  be  found 

2 

in  the  first  quarterly  report  . 

Consider  a  bent  plate  attached  to  a  circular  cylinder  as  shown 
in  Figure  1.  The  source,  a  monopole  antenna,  is  located  at  Q'.  The 
various  rays  contributing  to  the  radiation  pattern  are  illustrated 


1 


r 


Figure  1 .  A  bent  plate  attached  to  a  circular  cylinder. 

in  Figure  2.  Note  that  in  the  present  solution  higher-order  terms 
are  included  in  the  analysis.  The  source,  reflected,  and  diffracted 
terms  used  in  our  previous  solution^  are  shown  in  Figure  3.  The  super¬ 
imposed  pattern  using  these  three  GTD  terms  is  illustrated  in  Figure 
3d.  One  should  note  the  disconinuities  in  this  pattern  at  x  14°  and 
45°.  These  discontinuities  are  compensated  for  by  the  higher-order 
GTD  interaction  terms  illustrated  in  Figure  4.  As  before,  these  terms 
are  plotted  relative  to  the  same  signal  level  such  that  one  can  observe 
their  relative  significance.  It  is  clear  that  these  higher-order  terms 
can  be  significant  in  certain  sectors  of  the  pattern.  In  order  to 
demonstrate  how  the  complete  GTD  solution  creates  the  desired  total 
pattern,  various  combinations  of  terms  are  illustrted  in  Figure  5. 

The  geometrical  optics  solution  is  shown  in  Figure  5b  and,  as  before, 
several  discontinuities  exist  in  this  pattern.  The  total  superposition 
of  the  recently  introduced  higher-order  terms  is  shown  in  Figure  5c. 

It  is  very  interesting  that  the  patterns  illustrated  in  Figure  3d  and  5c 
superimpose  to  give  the  smooth  total  pattern  shown  in  Figure  5d.  As 
with  any  GTD  solution  for  a  complex  structure,  one  can  compute  even 
higher  order  terms.  However,  there  are  two  major  problems  with  adding 
such  terms  beyond  those  included  here:  1)  the  accuracy  of  the  GTD 
diffraction  solutions  become  questionable,  2)  the  numerical  solution 
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(c)  diffracted  field  (0)  (d)  S  +  R  +  D 

Figure  3.  Radiation  patterns  due  to  various  GTD  terms. 
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(a)  reflected/reflected 
fields  (R/R) 


(b)  reflected/diffracted 
field  (R/D) 


((.)  diffracted/ref  lected 
field  ( 0/R ) 


(d)  diffracted/diffracted 
field  (D/D) 


Figure  4.  Radiation  patterns  due  to  the  second 
order  interaction  GTD  terms. 


(a)  S  +  R 


(b)  geometric  optics 
solution. 

(S  +  R  +  RR) 


(c)  second  order  interaction  (d)  total  solution. 

GTD  terms. 

(R/R  +  R/D  +  D/R  +  D/0) 

Figure  5.  Radiation  patterns  due  to  various 
combinations  of  the  GTD  terms. 


becomes  very  inefficient.  Furthermore  based  on  cases  examined  to  date, 
it  appears  that  the  major  structures  found  on  aircraft  can  be  adequately 
solved  using  the  interaction  terms  considered  here.  In  order  to  illustrate 
the  various  GTD  mechanisms  they  are  all  shown  in  Figure  2  in  terms 
of  ray  paths..  The  previous  example  was  used  to  illustrate  the  various 
GTD  terms  included  in  the  present  analysis.  The  next  example  is  used 
to  demonstrate  the  validity  of  the  present  solution. 

The  geometry  illustrated  in  Figure  6  is  included  in  that  it  nds 
to  resemble  a  more  realistic  aircraft  configuration.  The  roll  p’  ie 
measured  and  calculated  patterns  are  shown  in  Figure  7.  Two  sets  of 
measured  and  calculated  patterns  are  illustrated  in  Figures  8  and  9 
for  68°  and  112°  conical  pattern  cuts.  In  each  case  both  the  Ep  and 
E^  components  of  the  near  zone  field  are  plotted  relative  to  the  same 
radiation  level.  In  addition  there  is  very  good  agreement  between 
calculated  and  measured  patterns.  The  above  patterns  are  presented 
to  illustrate  the  fact  that  our  numerical  solutions  can  solve  the  pro¬ 
posed  type  of  airborne  antenna  problem. 

A  second  aspect  of  our  analysis  of  airborne  antenna  patterns 
under  the  present  contract  is  to  investigate  ways  to  analyze  the  complete 
near  field  volumetric  pattern.  In  order  to  accomplish  this  goal,  the 
first  step  is  to  adopt  a  three  dimensional  model  for  the  aircraft  fuselage. 
The  model  must  provide  a  good  approximate  shape  for  the  fuselage  in 
the  antenna  studies  not  only  in  the  neighborhood  of  the  roll  and  ele¬ 
vation  planes  but  also  in  the  regions  close  to  the  nose  and  tail  of 
the  aircraft.  Furthermore,  the  model  must  be  efficient  and  of  a  form 
that  it  can  be  adopted  for  the  fuselage-wing  analysis,  i.e.,  finite 
flat  plates  can  be  easily  attached  to  allow  the  studies  of  fuselage¬ 
wing  interaction.  In  this  report,  a  surface  of  revolution  model  (namely, 
a  prolate  spheroid)  is  chosen  to  simulate  the  aircraft  fuselage.  Unlike 
the  previous  composite  cylinder  model,  the  prolate  spheroid  is  a  doubly- 
curved  surface  with  finite  principal  radii  of  curvature  and  is  apparently 
a  better  approximation  for  the  fuselage. 
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pattern  for  the  test  geometry  shown. 
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E4  .  68*  CONICAL  PATTERN 


ANGLE  (DEGREES) 


Figure  8.  Comparison  of  measured  and  calculated  near-field 
patterns  for  the  test  geometry.  Note  that  the 
engine  is  simulated  by  flat  plates. 
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RELATIVE  POWER  <dB> 


Figure  9.  Comparison  of  measured  and  calculated  near-field 
patterns  for  the  test  geometry.  Note  that  the 
engine  is  simulated  by  flat  plates. 


Once  the  model  is  decided  for  the  fuselage,  the  next  step  is 

to  investigate  the  radiation  from  antennas  on  this  doubly  curved  convex 

surface.  This  leads  to  the  studies  of  radiation  from  sources  on  general 

perfectly  conducting  convex  surfaces.  This  problem  was  treated  previously 

for  torsionless  surface  rays.  Subsequently  solutions  were  given  for 

3-5 

cylindrical  and  conical  surfaces  ,  where  torsional  surface  rays  exist. 
Under  the  present  contract,  the  previous  GTD  solutions^’^’^  were  extended 
to  obtain  a  uniform  high-frequency  solution  for  the  radiation  from 
apertures  and  monopoles  which  may  excite  torsional  surface  rays  on 
a  perfectly  conducting,  smooth,  convex  surface.  The  detailed  analysis 
of  this  work  was  presented  in  the  third  quarterly  report  ,  only  a  brief 
summary  is  given  here. 

Consider  a  source  located  on  a  perfectly  conducting  convex  surface 
at  Q',  as  shown  in  Figure  10.  The  high-frequency  solution  described 
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Fici.-rc  1C.  R:v?  emanating  from  a  source  on  a  convex  surface. 


here  employes  the  ray  coordinates  of  the  GTD.  Thus  in  the  shadow  region 
the  radiation  from  Q'  follows  a  surface  ray  to  Q  where  it  sheds  tangentially 
from  the  surface  to  the  field  point  P$;  whereas  in  the  lit  region,  the 
radiation  follows  the  incident  ray  of  geometrical  optics  to  in  the 
direction  of  the  unit  vector  s.  At  Q  the  orthogonal  unit  vectors  i 
and  n  (i.e.,  parallel  to  the  ray  and  normal  to  the  surface,  respectively) 
are  introduced  along  with  the  binormal  unit  vector  (b  =  t  x  n).  At  Q* 
these  same  unit  vectors  are  primed.  In  the  lit  region,  n  =  b  x  s  where 

A  /\ 

b'  =  b  is  perpendicular  to  the  plane  of  incidence. 

A  A 

The  high-frequency  electric  field  is  given  by  l  =  nEn  +  bE^  for 
points  away  from  the  convex  surface  in  both  the  shadow  and  lit  regions. 
Expressions  for  these  field  components  have  been  deduced  from  a  careful 
study  of  the  cylinder  and  sphere  conical  problems  in  which  higher  order 
terms  are  retained  in  the  asymptotic  solutions;  in  addition,  experimental 
results  for  sources  on  a  spheroid  were  helpful  in  the  generalization 
to  the  general  convex  surface. 

Consider  now  a  magnetic  current  moment  pm  tangent  to  the  surface 
at  Q'.  A  complete  description  of  this  solution  can  be  found  in  Reference 
6,  however,  a  simplified  form  is  presented  here. 


In  the  lit  region: 

■  -  f(p  . 

"n  4* 


'  =&  |(Pm-b,)(HS'+T^Fcos0i)  +  (pm.t,)TocoseiF]^1^  +  0(n£2)  (1) 


Eb  '  =${  (Pm’6’  >ToF+M'  )  (H^coseV) 


p" JkS  _p  o 

V- +  °( V’mJt  > 


(2) 
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In  the  shadow  region: 
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Note  that 


T0  -  Tpg(Q‘ ); 


ml  = 


kPg(q') 


1  1/3 


2(1+T2cosV) 


m  = 


kpg(q‘) 

~ T 


1  1/3 


F  =  S^-H^cose1 

l+T^cosV 


in  which  T  and  pg(Q')  are  the  surface  ray  torsion  and  radius  of  curvature 

at  Q',  respectively.  In  the  lit  region  the  various  terms  are  found 

by  projecting  the  incident  ray  onto  the  tangent  surface  at  Q* .  The 
0  0 

quantities  H,  H  ,  S  and  S  contain  Fock  type  functions  which  depend 
upon  the  surface  parameters  at  Q',  and  in  the  shadow  region,  upon  the 
surface  ray  trajectory.  The  angles  between  adjacent  surface  rays  at 
Q‘  and  Q  are  d\l>Q  and  d\|.’,  respectively  and  pc  is  the  caustic  distance 


In  the  monopole  case,  the  expressions  for  the  electric  field 
in  the  lit  region  are  given  as 

(5) 

(6) 

whereas,  in  the  shadow  region 
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+  0(m£  ),  and 


Eb '  ^  s"">V  HF +  0(mizi 


-jks 


-2i. 


+  0(m‘2),  and 


(7) 


(8) 


The  previous  expressions  reduce  to  the  geometrical  optics  field  in 
the  deep  lit  region  and  to  the  GTD  creeping  wave  field  in  the  deep 
shadow  region.  The  expressions  for  the  lit  and  shadow  regions  join 
smoothly  at  the  shadow  boundary.  As  expected,  they  reduce  to  the  asymp¬ 
totic  solutions  for  the  circular  cylinder  and  sphere  cases,  but  the 
higher  order  terms  in  "m"  must  be  retained  to  pass  smoothly  to  these 
two  limiting  cases.  As  the  radii  of  curvature  of  the  surface  become 
infinite,  TQ  =  0  and  Equations  (1),  (2),  (5)  and  (6)  simplify  to  the 
field  of  a  current  moment  on  a  ground  plane. 


Now  that  the  high  frequency  electric  field  is  given  in  terms 
of  the  surface  ray  coordinate  system  (namely,  a  ray  launches  from  the 
source  Q'  )  and  traverses  along  the  geodesic  path  to  the  shedding 
point  Q  ,  and  then  propagates  along  the  tangent  direction  at  Q  toward 
the  observation  direction  r(0r,<}>r);  the  major  task  remaining  is  to 
determine  the  unique  geodesic  path  for  a  given  radiation  direction. 

This  goal  is  accomplished  and  the  detailed  analysis  is  given  in  the 
second  quarterly  report7,  in  which  a  numerically  efficient  and  accurate 
scheme  has  been  developed  for  determining  the  geodesic  path  in  the 
case  of  an  antenna  radiating  from  a  general  convex  surface  of  revolution. 
The  surface  of  revolution  is  of  interest  in  that  it  provides  an  analytical 
model  for  the  aircraft  fuselage  structure.  A  computer  program  was  also 
developed  to  solve  the  governing  nonlinear  equations  using  the  secant 
(iteration)  method.  Some  numerical  results  in  terms  of  a  family  of 
geodesic  curves  are  presented  in  Reference  7  for  the  case  of  sphere, 
prolate  spheroid  and  cylinder. 
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Once  the  geodesic  path  is  determined,  one  is  ready  to  test  the 
newly  developed  solutions  (Equations  1-8)  by  applying  it  to  calculate 
the  radiation  from  slots  and  monopoles  on  perfectly  conducting  bodies. 

Some  numerical  results  are  obtained  for  the  cases  of  an  antenna  radiating 
from  circular  and  elliptic  cylinders,  cones  and  spheroids.  In  Reference  6, 
calculated  results  obtained  by  using  the  uniform  GTD  solutions  are 
compared  with  the  eigenfunction  solution  and/or  experimental  measurements 
and  the  agreement  is  excellent.  Some  examples  are  presented  here. 

The  patterns  of  a  radial  slot  on  a  cone  are  shown  in  Figure  11  which 


Figure  11.  Radiation  patterns  of  a  radial  slot  in  a  cone. 

compare  very  well  with  results  obtained  from  an  eigenfunction  solution. 

In  Figures  12  and  13,  the  patterns  of  a  circumferential  slot  and  monopole 
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Figure  12.  Radiation  patterns  of  a  circumferential  slot 
in  a  conducting  prolate  spheroid. 
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Figure  13.  Radiation  patterns  of  a  monopole  antenna 
on  a  conducting  prolate  spheroid. 

are  calculated  and  measured  in  the  plane  tangent  to  the  spheroid  at 
the  source  location.  The  prolate  spheroid  geometry  is  shown  in  Figure 
14.  Note  that  the  component  is  due  to  the  spheroid  surface;  it 
would  vanish  if  the  source  were  on  a  flat  ground  plane.  Finally,  a 
series  of  curves  for  the  b-component  of  the  electric  field  due  to  a 
circumferential  slot  on  the  spheroid  are  shown  in  Figure  15.  Notice 
that  the  calculated  and  measured  results  are  in  very  good  agreement 
for  the  cases  when  the  receiver  is  in  the  lit  region  (0=74°),  in  the 
shadow  boundary  plane  (3=90°)  and  in  the  shadow  region  (0=100°).  This 
confirms  the  validity  of  the  newly  developed  GTD  solution  for  an  antenna 
radiating  from  a  conducting  convex  body. 
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Figure  15.  Radiation  patterns  of  a  circumferential  slot 
in  a  conducting  prolate  spheroid. 
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II.  SIGNIFICANT  ACCOMPLISHMENTS 

A  series  of  accomplishments  have  been  made  under  the  present 
contract.  They  are  summarized  below: 

1.  A  paper  entitled  "Near  Field  Patterns  Analysis  for  Airborne 
Antennas,"  has  been  accepted  for  publication  by  the  IEEE 
Transactions  for  Antennas  and  Propagation. 

2.  An  oral  paper  was  presented  at  the  1978  International  IEEE/AP-S 
Symposium  at  University  of  Washington;  Seattle,  Washington. 

The  paper  was  presented  in  the  session  on  high  frequency 
diffraction  and  entitled  "A  Uniform  GTD  Solution  for  the 
Propagation  from  Sources  on  a  Perfectly  Conducting  Convex 
Surface," 

3.  A  paper  entitled  "A  Uniform  GTD  Solution  for  the  Radiation 
From  Sources  on  a  Convex  Surface"  is  under  preparation  and 
will  be  submitted  for  publication  in  the  IEEE  Transactions 
on  Antennas  and  Propagation. 
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